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I. INTRODUCTION 


A. MOTIVATION 


Binary computing systems have been used for about 30 years. The very-large-scale- 
integration (VLSI) technology has resulted in increasingly smaller circuits. However, 
the pace of reduction has showed because of two major problems, bus connection and 
pin limitation. In recent years, multiple-valued logic has been proposed to solve those 
problems [Ref. 1, 2]. Since the programmable logic array (PLA) is a basic tool for 
binary VLSI design, it is expected to be important in multiple-valued logic (MVL). This 
observation has inspired the work reported in this thesis. 

The truncated sum (TSUM) operation of multiple-valued logic is easily imple- 
mented in the charge-coupled devices (CCD) technology [Ref. 3, 4]. Several heuristics 
have been developed for the multiple- valued logic minimization involving TSUM and each 
claims some advantages in specific examples, but none of them is consistently better than 
the others [Ref. 5, 6, 7, 8, 9, 10]. Heuristic methods are interesting because of exact min- 
imization methods are extremely time-consuming. Until now, only Tirumalai and Butler 
[Ref. 9] have analyzed the performance of different MVL minimization algorithms. Their 


study is limited to two-variable four-valued functions. Three open questions exist: 


1. In the case of two-variable four-valued functions, can we design an algorithm which 
is better than previously proposed algorithms? 


2. What is the algorithm performance if we increase the number of variables (greater 
than two)? 


3. What are the computation time requirements for the various algorithms? 


With the computer software developed at Naval Postgraduate School called HAM- 


LET (Heuristic Analyzer for Multiple-valued Logic Expression Translation) [Ref. 10], 


users can easily investigate their own heuristics. For example, the neighborhood de- 
coupling (ND) algorithm described in this thesis is built as one independent option of 
HAMLET. This thesis reexamines the algorithms used in HAMLET and analyzes the 
new method, the ND algorithm. It also investigates each algorithm under four different 


function settings and summarizes their performance and computation time. 


B. BACKGROUND 

The heuristic approaches in HAMLET can all be classified as greedy algorithms. 
The HAMLET execution procedure of these algorithms is abstracted as follows. Formal 
definitions will be covered in the next chapter. Let f be a multiple-valued function and 


let a, a minterm, be an assignment of values to the variables of f such that f # 0. 


je 2K OK KKK KKK kK KK KK KK IK KK KK KK CIE 2k I IC oo KC KC IK I I KC IK IC 2K IK IK 2K KC KC KK OK KK 2k IK 2 OK 2K OK OE KK 2k 


Input: let the Mf be the set of minterms of a function f; 


Output: the minimized sum of product, S, of the original function; 


2K KOK OK OK OK KK OK KK KK KK KK KK KC KK IE IC IC KE I 2k KI OK IK I IK IK IK IK I OE KKK IK 2K KK 2K 2K 2K OK KKK KK 2K OK ~h 


{ 


S<— @. 
While (Mf 4 ¢) do { 


pick one minterm a from M; 

find an implicant J, which covers a; 
S—1I,US; 

subtract [, from f; 


} 


TABLE 1.1: Summary of Three Heuristic Algorithms 


Heuristic Algorithm Choice of Minterm Choice of Implicant 


Pomper and Armstrong [Ref. 5] Random Drives Most Minterms to 
pera) |_| Dor dover 
Besslich [Ref. 6] Smallest Weight | Drives Most Minterms to 
0) 
Dueck and Miller [Ref. 7, 8] Largest IF Largest BCR 
88 Loton tates) | 


Table 1.1 shows three previously proposed algorithms. Each algorithm differs from 













the others in the manner of picking the minterms (qa) and finding the implicants (J,). 
For example, in the above recursive search procedure, the input function expression f is 
evaluated at minterm a. Next an implicant J, is chosen which covers a. Then, implicant 
I, is added to output solution set S, and is subtracted from function f. 

Pomper and Armstrong introduced a heuristic algorithm that picks a randomly (as 
long as a is in the set of minterm A) and finds J, (as long as J, covers a) which drives the 
most minterms to 0 or don’t-care when I, subtracted from function f [Ref. 5]. In 1986, 
Besslich presented a direct cover algorithm, according to weight transformations. The 
Besslich algorithm picks a with the smallest weight (most isolated minterm) and finds I, 
which has a lowest cost per minterm covered (i.e., which drives the most minterms to 0 
or don’t-care) [Ref. 6]. In 1988, Dueck and Miller presented another algorithm that picks 
a from M if a has the highest isolated factor (IF) and then finds the J, which directly 
covers @ such that the break count reduction (BCR) is maximum [Ref. 7, 8]. 

Fach algorithm has some advantages in specific examples [Ref. 9]. The ND algo- 
rithm is characterized by adopting the advantage from each algorithm and fully utilizing 
the properties of the truncated sum. The ND algorithm is an improvement to the Dueck 
and Miller algorithm with revised decision rules for making selections of minterms and 


implicants. 


C. THESIS OUTLINE 

A summary of MVL definitions for truncated sum minimization are introduced in 
Chapter II. The notations and definitions of Chapter II also help us in explaining the 
algorithms appearing in subsequent chapters. The neighborhood decoupling algorithm 
is presented in Chapter III. Chapter IV and V discuss the performance and computa- 
tion times of the neighborhood decoupling algorithm with Pomper and Armstrong’s and 


Dueck and Miller’s algorithm for different function parameters. 


Il. NOTATIONS AND DEFINITIONS 


Although multiple-valued logic minimization for truncated sum has been discussed 
in resent years [Ref. 7, 8, 9, 10, 11], a distinct and complete formal definition for truncated 
sum MVL minimization is not available. In view of this, we studied the truncated sum 
MVL minimization, then summarized and give several definitions that will be presented 
in Section A. It will also help us to explain the algorithm given in the next chapter. We 
will report our observations on truncated sum properties in Section B, and the definitions 


for ND algorithm will be discussed in Section C. 


A. DEFINITIONS FOR TRUNCATED SUM 


Definition 1: 


Let ¥ = { 21,22,.-.,2n } be a set of n input variables where z; takes on values 


from R = {0,1,...,r —1}. An n-variable r-valued function f is a mapping 
f:R" — Ru {r}. [Ref. 9] 


Here, r is a don’t — care value; it can be chosen freely from any of the logic values, 


Definition 2: MIN 


The MIN [Ref. 9] function is denoted as f(x,,x2) = 2,r2 which evaluates to the 
minimum value of its arguments. For example, if R = {0,1,2,3}, then f(1,2) = 1 and 


f(0,3) = 0. A minterm is an assignment of values to 21, 72,...,2, such that f(x) ¥ 0. 


Definiton 3: Literal 


The literal operation of a variable z is defined as: 


ab _ r—1 ifa < MY = b 
‘ = 0 otherwise. [Ref. 11] Ca) 
Definition 4: Truncated Sum (TSUM) 
The truncated sum (TSUM) operation is defined as: 
TSU M(a@y, 22) = 21 + 22 = aria,  2o,r — 1). [Rey 11 (2.2) 


The two + signs in this expression are different. The leftmost denotes the TSUM 
operation, while the rightmost denotes ordinary add two of logic values which are viewed 
as integer. The TSUM obeys the associative and commutative rules. For example, if R 
= {0,1,2,3}, then TSUM(1,2)=s andepstr a7 j—2) 

These definitions are inspired by the fact that CCD implementation supports TSUM 
naturally (Ref. 3, 4]. 


Example 1: 


For example, 'z} is a literal and takes value of 3 when 1 < x, < 3. However, 
function 2 ‘x? takes a value of 2 based on to the definition of MIJN. Similarly, the 


product term 2 323 °z9 is a function that is 2, when z, = 3 and r2 = 0. 


Definition 5: Product Term 


A product term p is the MIN of one nonzero constant c € 7, and one or more 


literal functions. In general, a product term is defined as: 


ie epee in», Jn tk S Jk 
9 = eon ree ben i Sree (233) 


The constant or coefficient c, in a product term, effectively scales the term. For 
each variable z;, we say the window size of the literal **a;2* is 7, —1, +1. We use the 


terms product term and implicant interchangeably in this thesis. 


Definiton 6: Sum-of-Products Expression 


A sum-of-products expression is py + po +...+ pn for some integer N, where p; is 


a product term. 
Definition 7: Minterm 


Pemminterm a is of the form ez,’ 725° ~a2,° ...°"a-” where a; € R and constant 
c€ R— {0}. We say the coordinate of a is < a),a2,...d, >. We denote the value of 
minterm a, g(a), as the nonzero constant c. 

A product term p = c "2,7! ? x27? ...'"2,2" can be decomposed into []7_, (jx —7e +1) 
minterms. We say p generates those minterms. Given a product term p, the set of 
minterms generated from p is denoted by A/S,. If the number of elements in J/S,, is 
greater than that in M[S,,, we say p; covers a larger area than p2. Given a function f, 


the set of minterms generated from its product terms is denoted by MS. 
Definition 8: Saturated Minterms (SAT) 


Given a minterm a generated from the original function to be minimized, if g(a) = 
r—1l1, then a is a saturated minterm. Let SAT be the set of all saturated minterms of a 


function. 
Example 2: 
If the input function to be minimized is expressed as follows, 
f = 3 %rp'2, + 2 bapa, +3 bay *a2 + 2 %ay*ay + 1 Sap ory + 1 zy 2g, 


the MS; can be represented as 11 minterms in Figure 2.1. We mark a saturated minterm 


with a dot in the figure. 





Figure 2.1: Map for Example 2, 3, 4; Step 1 of Table 3.2 


Lemma 1 Given a minterm a the mazimum number of implicants which covers a 1s 


Chg). 


Proof: Consider a variable (axis) 2; of a. Any implicant (/,) that covers a may have a 
range or “window size” w, such that 1 < w <r. With a window size w, we may have w 
implicants that covers a. That is, the minterm a@ can be at any position of the window. 
For a given axis, counting all window sizes, we have 1+2+3+...+r= arth) possible 


implicants that cover a. Over the entire n-dimensional space, we have (cee) = "Ore 


[Ref. 12]. 


B. THE PROPERTIES OF TRUNCATED SUM 


There are two important properties of the truncated sum which are useful later in 


developing the ND algorithm. 


1. Saturated minterm can be generated by TSUM operation. 

The truncated sum of two or more minterms may produce a saturated 
minterm. By definition 4, the truncated sum of any minterm and a sat- 
urated minterm remains a saturated minterm. In other words, given two 
minterms a, 8 and minterm y = TSUM(a, (3). If value of 7 is r —1, i-e., 


y is a saturated minterm then for any other minterm 6, y+ 6 = 4. 


As an example, in a 2-variable 4-valued function, three minterms add 


in one position. 
2 
2 Uyl 2g? 49 Vy) 2p2 4] Vyl 2p? = 3 gl 2p? 4 1 Vy) 2g? = 3 ly) 2y? 


The first two terms form a saturated minterm, and this saturated minterm 


absorbs the third minterm. 


2. Don’t care minterm can be produced by saturated minterm. 

In the minimization procedure, we may update a minterm a to a’ 
by subtracting minterm y (a’ = a — y), where 7 is the value of selected 
implicant. If a € SAT, in succession of updates, the value of a’ may reach 
the value 0. In that case, the algorithm will reset that minterm coordinate 


to don't-care, 1.e., value r. 


C. DEFINITIONS USED IN ND ALGORITHM 


Definition 9: Direct Neighbors 


Let a and @ be minterms with coordinates < a,,a2,...a, > and < bh), b2,...6, > 
respectively. If for all 2 we have a; = 6; except one position 7 such that |a; — 6;| = 1 we 
say that a and @ are direct neighbors. Given a minterm a, we use N(q@) to denote the 


set of its direct neighbors. 
Observation 1: The mazimum number of direct neighbors of a given minterm is 2n. 
Definiton 10: Directional Neighbors 


Two minterms a and £ are directional neighbors in the direction z;, if a; = 6; for 
all 2 € [1,n] such that 7 # 7. When 6; > a; we say that § is in the positive direction of 


a, while when b; < a; we say that £ is in the negative direction of a. 





Figure 2.2: Example for Connected Minterms 


Observation 2: /f 8 is a direct neighbor of a then B is a directional neighbor of a in 


the direction of x; for some i € {1, nl. 
Definition 11: Connected Minterms 


This is a recursive definition. Given a minterm a and a minterm £, then we say 


is a connected minterm of a, if 


1. B is a direct neighbor of a and either g(3) < g(a) or a € SAT. 


2. 8 is a directional neighbor of a in direction z; and (’s direct neighbor is 
connected to a and either g(8) < g(a) ora € SAT. 


20.0 3731 1 Bec 


For example, in Figure) 2:2 tminterinmess 0 le v3 and 2027 


(pointed by arrow) are connected minterms of 2 ?z?/z} (the minterm with @ sign). 
Definition 12: Connected Minterm Count 


CMC, is the connected minterm count of minterm a. It is the number of minterms 


that are connected to minterm a. 
Definition 13: Expandable Directional Count 


EDC, is the expandable directional count of minterm a. It is the number of 
directions (both positive and negative for each z;) in which @ has one or more connected 


minterms. 


10 





Figure 2.3: Map for Example 3, 4; Step 2 of Table 3.2 


Observation 3: 0 < EDC, < 2n. 
Definition 14: Clustering Factor 


The clustering factor relative to a minterm a is defined as 
CF, = EDC.(r ~1) + CMC,. (2.4) 


This is a measure of the weight of all connected minterms relative to a. The (r—1) 


factor is the range, or maximum possible number of minterms, in a direction 2;. 


Example 3: 


In Figure 2.1 the minterm 1 °r°*z? (the minterm with * sign) has no connected 


minterms nor expandable directional neighbors, i.e., its CA{C and EDC values are 0. 
Figure 2.3 shows that a circled implicant 1 °r}!z3 was subtracted from Figure 2.1. We 
mark a minterm with a dot in figure because it was a saturated minterm in original 
function map (see Definition 8 and Figure 2.1). The minterm a = 2'zr} °r$ (the minterm 


with @ sign) is one of the ten minterms and CMC, = 3, EDC, = 2. The clustering 


factors of all minterms in Figure 2.3 are listed in Table 2.1. 


Wal 


wee 2. af Es oiels eu Sa rsehents in Fig ue jay a 
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[Minterm | 372° %eqh [Oa Pay? [OPP ay? | IIa ye Saya 
[or | 6 [| @ | ¢ [| 0 | 0 










Vz 


Ill. NEIGHBORHOOD DECOUPLING 
ALGORITHM 


The neighborhood decoupling (ND) algorithm [see Appendix A for the C program 
listing] is an improvement to Dueck and Miller’s method [Ref. 8]. Like the algorithms 
mentioned in Chapter I Section B and throughout, the ND algorithm has two compu- 
tational phases: select a minterm and select an implicant. Firstly, the most isolated 
minterm is chosen by using the algorithm M to be described in Section A. The most iso- 
lated minterms in general are different from Dueck and Miller’s method due to different 
decision rules. Secondly, from all implicants which cover the most isolated minterm ND 
algorithm chooses the one that is not strongly “coupled” with its neighbors. In other 
words, it is the most isolated implicant. This decoupling process is based on the ob- 
servation that if we choose that specific implicant then we may minimize the negative 
impact for future minterm selections as well as implicant selections. The explanation of 
this idea is described in Section B. In the algorithm below, f denotes the function to be 


minimized. 


{ 
SS — ¢; /* SS = Solution Set */ 
WS = MS; = {ala is generated by the function f; if a € SAT then mark its 
coordinate }. 
While WS £ ¢ do { 
1. Use algorithm M (see Section A) to select a minterm a from the WS. 


2. Use algorithm N (see Section B) to select an implicant J, that covers a. 


See) eee Oe Pe 
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4. VB EI, do { 
compute g(8) — g(f) — g(a). 
subtracted J, from WS. 
if GB is originally marked and g(?) = 0 then g(Z) — r. 


/* don’t-care terms */ 


A. ALGORITHM M: MINTERM SELECTION 
Based on Definition 15, all minterm’s clustering factor can be computed. The ND 
algorithm computes the clustering factor in the order of coordinates (x;), 1.e., row-column 
order. For example, the minterm 2 *z? °z9 is evaluated earlier than the minterm 2 7x7 723. 
The algorithm M is described as follows: 
1. Compute the corresponding C’Fy, for all a, Ee WS. 
2. Select the minterm a that has the smallest clustering factor (excluding the minterm 
a € SAT as well as greater than 2 x (r —1)). If there is a tie, the first one that 
gets evaluated is chosen. 
B. ALGORITHM N: NEIGHBORHOOD RELATIVE COUNT 
The purpose of algorithm N is to choose the most “isolated” implicant (/,) and 
update the working set WS. It computes the neighborhood relative count (NV RC) for all 
implicants that cover the minterm a. The implicant with the smallest N RC is chosen. In 
other words, NRC is a measure of the coupling strength of an implicant with its neigh- 
bors. To select an implicant which is equivalent to breaking the coupling between that 
implicant with its neighbors, the candidate implicant should have the smallest coupling 


strength with its neighbors. Therefore, the ND algorithm tends to choose the most “iso- 


lated” implicant. If there is a tie in selecting the [,, ND algorithm chooses the one which 
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covers the largest area. The computation of NRC for a given implicant is described as 


follows: 
1. Initialize the NRC to zero. 


2. Check all neighboring minterms of the implicant and increment or decrement its 
NRC according to the following (intuitively stated) rule, which is, if the coupling 
strength between covered and uncovered area is weak (good for further decoupling), 
algorithm N decreases N RC’, otherwise increases N RC. 


i AE KK KK OK 2k KK KK OK KK 2K 2K KK KK IK IK IK I IK OK og IK KK OK IK IK OG IG IG IK OI IK I IK I I I I IK OIC IK 2 I I I OK IK IK IK I IK IKK IK I IK IK IK OK IK Ok 


a: the chosen minterm from algorithm M. 
M: the set of minterms which was covered (generated) by the chosen impli- 
cant (J,). 


N(8): the set of direct neighbors of minterm £. 


2K OK IK 6 IK KK KK KK I KK IK IG OK IK IK IK IK IK IK I IC IK IK IK IK IG I I IK og IK IR IC IG IK IK IG OIC AG IK IK IC IK OIC OE IK I I I OK fg IK IK I oS IK OI OK OI IK OK OK ote a 


{ 
NRC «< 0; 
V8 eM and B 4a do { 
if(g(B) — g(a) < 0) then NRC — NRC — 2: 
} 
V3 € M and Vy € N(8) do { 
if (y g M and y £0 and (7 ¢ SAT or 8 ¢ SAT)) then { 
if (9(8) — g(a) > g(y) ) then { 
if (y € SAT ) then NRC — NRC - 1; 
else NRC — NRC +2; 


if (g(8) — g(a) < g(y) ) then { 
if (g(8) = 9(7) ) then NRC — NRC + 2; 
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if (y € SAT and g(y7) — g(8) <0 ) then 
NRC — NRC +2; 

else { 

if (9(8) > g(a) and 9(8) ¥ 9(7) ) then { 
if( 8 € SAT) then NRC — NRC - 1, 

else NRC — NRC +2; 
} /* end if */ 
} /* end else */ 
Le endaiaay 


if (9(3) — g(a) = g(7) ) then { 
if (g(y) > 0 or B € SAT) then 
NRC — NRC -1; 
else NRC — NRC — 2; 
} 
} /*end if */ 
} 
if (M = {a}) then { 
if (a € SAT) then NRC < 2; 
else if (NRC <0) then NRC <— 1; 


} 
else NRC — NRC +2; 
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Figure 3.1: Map for Example 4; Step 3 of Table 3.2 


Example 4: 


It is instructive to examine an example of the application of the ND algorithm. 
We use a two-variable four-valued function as an example. The input function to be 


minimized is expressed as: 
f= 3S tal 4 214279 4 3 Vy 293 4 2 2p? 293 4 1 2923q3 4 1 9g 2y?, 


The working set, WS, is initialized to MS; and is represented in Figure 2.1. The clus- 
tering factors of all minterms in WS are calculated (see example 3 for computation). 
The smallest CF that first comes from minterm 1 °r??z2 ; therefore, algorithm M will 
select a = 1 °r9*r2. The ND algorithm computes the NRC, for each implicant J which 
covers @ using algorithm N. Since implicant 1 °x}!z? has the smallest NRC(—2), the 
ND algorithm selects it as the first implicant in the solution set (SS). Table 3.2, together 
with Figures 2.3, 3.1, and 3.2 shows the steps of choosing successive implicants. The 
* sign in each figure indicates the most isolated minterm while a circled implicant is 
the most isolated implicant. Suppose we have chosen two implicants from the function 


shown in Figure 2.1 the resulting function is shown in Figure 3.1. The minterm 1 7x? %z3 


is selected, since it has the smallest C'F'(4). 


Uf 





Figure 3.2: Map for Example 4; Step 4 of Table 3.2 


TABLE 3.1: NRCs for Minterm 2 °r)'r) in Figure 3.2 
implicant [372,92 ary 
DwRO [= fo 





There are two implicants that cover the minterm 1 *r?%z73, and their NAC values 
are 2 and —2. The implicant 1 'r?%z3} is chosen since it has the smallest NRC(—2). 
Having updated the working set and added 1 !z?%z3 to the solution set, we have the 
new map in Figure 3.2. From Figure 3.2, the NRC values for the minterm 2 °r?'!z} are 


available in Table 3.1. Finally, the working set should contain value 0 (empty square) or 


4 (don’t-care) as shown in Figure 3.3. 





Figure 3.3: Final Working Set 
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TABLE 3.2: Steps of ND algorithm 


Step | Minimum | Minterm a | Minimum | Candidate Implicant | Reference 
a — 


2 ee ais oe oF 


a a 
ear 






The final minimized result, g, 1s expressed as: 


g=2%z 1 3 +1 ‘2% 20) e+". t,t +3 ows 31 Lp. 
2 2, 


Is 


IV. COMPARISON RESULTS 


In this thesis all testing results were obtained by running sample functions on the 
VAX 11/785 and ISI workstations. A large number of sample functions were randomly 
generated (63,500). We applied each algorithm to minimize these sample functions in 
a way similar to the method used by Tirumalai and Butler [Ref. 9], and Yurchak and 
Butler [Ref. 10]. However, Tirumalai and Butler generated the input functions with 
a fixed number of minterms. We generated the input function with a fixed number of 
product terms. This thesis investigated three algorithms: (1) Pomper and Armstrong 
[Ref. 5], (2) Dueck and Miller [Ref. 7, 8, 10], and (3) Neighborhood Decoupling under 
various settings: 

1. Two-variable four-valued with 3 to 16 input product terms. 
. Two-variable five-valued with 3 to 25 input product terms. 


. Three-variable four-valued with 3 to 45 input product terms. 


=> OW bd 


. Four-variable four-valued with 3 to 35 input product terms. 


For two-variable four-valued functions, we generated 1000 sample functions for a 
given number of input product terms from 3 to 16, 1.e., we generated 14,000 functions 
to test the different algorithms. For the other three, we generated 500 functions for each 
given number of input product terms, 1.e., for each set we generated 11,500, 21,500, and 
16,500 functions respectively. 

In Section B, functions with a different number of product terms are tested and the 
computation times are recorded. For a given number of product terms, the computation 
time from two to four variable functions grows exponentally. This explains why we did 


not produce more than 1000 functions for each case. 
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A. PERFORMANCE COMPARISON 


There are two ways to measure the performance. First, the average number of 
output product terms will show us the advantage of each algorithm. Second, the ratio 
of the performance results will show us the performance increases or decreases for each 
algorithm under different conditions. 

1. AVERAGE NUMBER OF OUTPUT PRODUCT TERMS 

For each set, we compute the average number of output product terms (see Ap- 
pendix B). From these data points, we can plot a curve, indicating the average number 
of output product terms as a function of input product terms. The plots for the various 
sets are shown in Figures 4.1, 4.2, 4.3, and 4.4. In these bell-shaped figures we observed 
three important features: 

1. The differences in the number of average output product terms among those algo- 
rithms are insignificant, when the input number of product terms are small. When 
the number of input product terms is small, the function is simple and easy to 
minimized; therefore, the number of output product terms for each algorithm are 
almost same. 

2. The highest point of each bell-shaped curve is the highest average number of output 
product terms for each algorithm. The number of input product terms that has the 
highest average number of output product terms constitutes the most complicated 
functions. 

3. When the number of input product terms get larger, the curve goes down gradually. 
This is chiefly due to the fact that the more input product terms in a function the 
higher tendency of generating saturated minterms. In most cases, a single implicant 
can cover a Cluster of saturated minterms. 


From Figures 4.1, 4.2, 4.3, and 4.4, the neighborhood decoupling algorithm obvi- 


ously outperforms the other two algorithms. 


a) 


Average Number of Output Product Terms 


Average Number of Output Product Terms 
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Dueck & Miller 
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Figure 4.1: Two-Variable Four- Valued Average Product Term 
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Figure 4.2: Two-Variable Five-Valued Average Product Term 
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Figure 4.3: Three- Variable Four-Valued Average Product Term 
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Figure 4.4: Four- Variable Four-Valued Average Product Term 
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TABLE 4.1: Test Results of 14000 2-Variable 4-Valued Sample Functions 


Pomper and Armstrong | Dueck and Miller | Neighborhood Decoupling 


best 94 212 336 
better 479 2005 2308 
equal 11002 11002 


2. PERFORMANCE RATIOS 














We consider another measure, the performance ratio, to demonstrate the perfor- 
mance of each algorithm. First, for each algorithm, we count the number of minimized 
functions that specify which algorithm is the “best” of the three, i.e., the number of in- 
stances that a specific algorithm uses the minimum number of implicants (output product 
terms). If two algorithms use an equal number of implicants and less than the other one, 
we determine that they are “better” than the third one. When all three algorithms use an 
equal number of implicants to minimize a function, we determine that they are “equal”. 


The performance ratio is defined as 


le [Noes a aNiberien =r Nese) as ae (45 


Where Neest, Noetter, and Neguat are the number of instances that specific algorithm per- 
forms “best”, “better”, and “equal” respectively. 

The total number of functions tested, Nota: in our case is 14,000, 11,500, 21,500, 
or 16,500. Tables 4.1, 4.2, 4.3, and 4.4 show the performance ratios for each setting. 
For example, in Table 4.1, with 14,000 functions tested, we counted the cases which 
Neighborhood Decoupling algorithm performs no worse than the others as 13,646. That 
is, Wnp = 13646/14000 = 0.9747. Tables 4.1, 4.2, 4.3, and 4.4 show that the performance 
is degraded when n or r is increased. However, the neighborhood decoupling algorithm 


obviously outperforms the others. 
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TABLE 4.2: Test Results of 11500 2-Variable 5-Valued Sample Functions 
best 114 396 
better 
equal 


8026 9499 10907 
0.6979 0.8260 0.9484 













TABLE 4.3: Test Results of 21500 3-Variable 4-Valued Sample Functions 


Pomper and Armstrong | Dueck and Miller | Neighborhood Decoupling 


948 1595 
2791 4307 
(fol 1871 


11610 13773 18367 
0.5400 0.6406 










TABLE 4.4: Test Results of 16500 4-Variable 4-Valued Sample Functions 


Pomper and Armstrong | Dueck and Miller | Neighborhood Decoupling 















best 1604 S0nt 
better 1420 3106 
equal 2884 2884 


0.3581 0.5495 0.6834 
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B. TIMING COMPARISON 

HAMLET CAD (Computer Aided Design) tool has an option that can measure 
program’s computation time. In this section, the algorithm timing comparison counts the 
average computation time for each input product term minimization. These measures 
are collected on the ISI workstations at NPS. The theoretical worst case computional 
complexity of each algorithm is O(r?") (see Chapter II). The decision rules for selecting 
minterm a and implicant J, in each algorithm are different. In general, an algorithm 
which needs more or complex decision rules in selecting a and J, should take longer 
computation time. There may be a trade-off between performance and computing time. 
In each computation cycle (select a minterm and select an implicant), each algorithm 
generates an output product term. The ND algorithm uses more complex decision rules 
than Dueck and Miller’s. It may appear that ND algorithm should use more computation 
times. However, this is not true, since ND algorithm in general has less computation cycle 
than other algorithms so that it stop computation earlier than other algorithms. From 
the analysis of testing case results (see Appendix C), the average computation time of 
the ND algorithm is no worse than Dueck and Miller’s for more than 4 second, sometimes 


even better. 
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V. DISCUSSIONS AND CONCLUSIONS 


A. DISCUSSIONS 


Notice that in the comparisons we did not include the Bessilich’s algorithm and the 
absolute minimum solutions. We now justify the reasons. Besslich’s algorithm was not 
included in HAMLET because it requires a truth table to operate efficiently. We did not 
have room to store large truth tables. Also Besslich’s algorithm 1s likely to require much 
time [Ref. 10]. In addition, we know that the Dueck and Miller is a satisfactory [Ref. 9, 
10] heuristic which is reasonable approximation to the absolute minimum solutions for 
two-variable four-valued functions. 

The performance results show that the ND algorithm outperforms the other two 
algorithms. Although Pomper and Armstrong’s algorithm is faster than the other two 
algorithms, about one third of the computation time, the number of average output 
product terms and performance ratios are worse than the other two algorithm’s. From 
the timing comparison results, the ND algorithm runs faster than Dueck and Miller’s in 
larger variable setting conditions (more than two variable). This time efficiency is due 
to the decision rules employed in the ND algorithm which take advantage of the special 
property of truncated sum operations. In other words, the input product terms have a 
tendency to produce saturated minterms. In the decoupling process (algorithm N), a 
minterm in SAT will always qualify to combine with its neighbors to form an implicant. 
Although the ND algorithm conceptually is similar to Dueck and Miller’s algorithm, the 
ND algorithm uses saturated minterms in an more effective way. For example, when 
we update (deduct) saturated minterm from expression, the minterm will be updated to 


a “don’t-care minterm” (see Chapter III). Like binary logic minimization, a “don’t-care 


Zz 


minterm” [Ref. 13] can simplify the minimization process. 

Recall that in algorithm N, we compute the NAC values for a given minterm @ by 
examining the relationships of J, and its immediate neighbors, 1.e., one step look-ahead. 
It is natural to believe that with more steps look-ahead we might make a better choice 
of the implicant and therefore provide a better solution. The exponential growth of the 


number of possible implicants restricts the practical use of k-lookahead for & # 1. 


B. CONCLUSION 

The truncated sum MVL minimization can be done by the neighborhood decou- 
pling algorithm which selects the most isolated minterms as well as implicants. In the 
development of ND algorithm, we have the following observations: (1) truncated sum 
operations may produce saturated minterms, and (2) a saturated minterm will reduce to 
a don t-care minterm in the minimization process. The decision rules of the ND algorithm 
take full advantage of these observations. The ND algorithm outperforms most heuristic 
methods and does not lose its run time efficiency because the algorithm finds the solution 
and stops earlier than others. As the number of variable of the input function become 
larger, 1.e., greater than three, the ability to minimize the function in reasonable time is 


an important and challenging research area. 
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APPENDIX A: PROGRAM LISTINGS 


Paend.c 


- This module implements the Neighborhood Decoupling heuristic for 
MVL minimization. 


- This program can be added to HAMLET [Ref. 10]. 


- HAMLET is written by the instructor J. Yurchak in the Department 
of Computer Scinece at Naval Postgratuate School. 


* | 
#Hinclude "defs.h" 
Seacic int 


better_found; 


feaieighborhood Decoupling functions ~-~-<-~-<-<-<--<<-<<-<-<<<<<<<<<<-<<------- * / 


N_D() 
[rn ee ee ee ee = 
eeunction: 
- Perform the Neighborhood Decoupling heuristic on the input expression 
:algorithm: 
Start with a working copy E_work of the original function E_orig; 
Initialize a final function E_final; 
While (there are still minterms to pick) f 
Pick a minterm X from E_work; 
Pick the best implicant I for X; 
Add I to E_final; 
Subtract I from E_work; 
} 
:globals: 
E_orig 
e_flag 
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m_flag 
q_flag 
G_flag 
FO_ratio 
:Side_effects: 
STAT 
HEUR 
E_work 
E_final[(] 
(Gar Oda by. 
main() 
-Cdlice 
dealloc_expr() 
dup_expr () 
print_terms () 
print_map() 
mim() 
pick_implicant() 
subtract_implicant() 
print_source() 


{ 

int num_impl = 0; 
191 *X ; 

Implicant *I; 

float Bac le; 


if (E_final(ND].I != NULL) 
dealloc_expr(&E_final [ND]) ; 


#Hifdef ANALYZER 
STAT = &WY_stat; 
#Hendif 


HEUR = ND; 

dup_expr(&E_work, &E_orig) ; 
E_final{HEUR] .nterm = 0; 

E_final (HEUR] .radix = E_orig.radix; 
E_final(HEUR] .nvar = E_orig.nvar; 
E_final{HEUR].I = NULL; 


#ifdef ANALYZER 
if (e_flag) 
print_terms(&E_orig) ; 
if (m_flag) { 
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Print. Grmig map (nda) va), 
print_map() ; 
} 
#endif 


better_found = Sc_flag; 
resource_used(START) ; 


for C35) 
aoe -eminGleework))) == NULL) { 
if (num_impl < E_orig.nterm) 
better_found = 1; 
break; 
4 
I = pick_implicant(X); 
num_impl++; 
subtract_implicant (I) ; 


#ifdef ANALYZER 
eiGuet lag) 
print_implicant(X,I); 
if (m_flag) 
print_map(); 


#endif 
ie (Smit Lap) 
if (num_impl >= E_orig.nterm) 
break; 
i 
ii 


resource_used(STOP) ; 


if (!better_found) { 
num_impl = E_orig.nterm; 
dup_expr(&(E_final[ND]) ,&E_orig) ; 


#ifdef ANALYZER 
if (!q_flag @& !G_flag) { 
if (!better_found) 
printf(" %-4d nda: %4d /%-4d %61d.%2.21d\n\n", 
expr_seq,num_impl,num_impl,secs_used() ,tsecs_used()) ; 
else 
printf(" %-4d nda: %4d /%-4d %6d.%2.21d\n\n", 


33 


expr_seq,num_impl,E_orig.nterm, 
secs_used() ,tsecs_used()); 


Ip 
#Hendif 


dealloc_expr(&E_work) ; 


Static vin *xmim(E) 
Expression *E; 


[R pn one ne 


7 Lune ton: 
- Compute the clustering factor and find the Most Isolated Minterm 
in the expression pointed to by E, and return its coordinates as 
a vector. 
- Local to nd.c 
:globals: 
radix 
nvar 
:side_effects: 
STAT 
-called by: 
N_D() 
calls: 
next_coord() 
eval() 
vcopy () 
:returns: 
- A vector of integers representing the coordinate of the most 
isolated minterm, or NULL if no more minterms. 
- The value at that location is also returned as the last integer 
in the vector. 


register chy These 

int cur_val = E->radix, 
cur_CF = MAX_INT, 
X_orig [MAX_VAR+2] , 
Ro = radreeat, 
Not_all = QO, 
All_trun = OQ, 
TRUN = 2*R_1, 
last = QO, 
expanded, 
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value. oie 
vali(2}, 
Viel 2 lize 
Er 
ea, 
dea, 
term; 
int *X,*next_coord(); 
Seacic int 
coord[MAX_VAR+2] , 
save_coord [MAX_VAR+2] ; 


#ifdef ANALYZER 
STAT->calls_mim+t+; 
#endif 


for (term=0; term < E_orig.nterm; term++) { 
k = 1; 
while ((X=next_coord(coord,&(E->I[term]),k)) '= NULL) { 
vcopy(value,eval(E,X)); 
if (value[EVAL] && value[EVAL] < radix) { 
if (!value[HLV]) 


Not_all = 1; 
if batrun) 4 
Gi = 0, 
dea = 0; 
ea = QO; 
for (j=0; j < nvar; j++) X_orig([j] = x(jl; 
/* for each variable (direction)... */ 


for (j=0; >) <nvar; jr+ ) 
expanded = 0; 
/* If not on a left hand edge, move left */ 
while (X[j] > 0) f 
Xijlea 
vcopy(vali,eval(E,X)); 
if (vali[EVAL]) { 
expanded = 1; 
eat+; 
} 
else 
break ; 
} 
X{j] = X_orig(ljl]; 
if (expanded) { 
expanded = 0; 
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dea++; 
} 
/* if we didn’t start on a right hand edge, move right */ 
while (X[j] < R_1) { 
Mig dite 
vcopy(val2,eval(E,X)); 
if (val2[EVAL]) { 
expanded = 1; 
eatt; 
} 
else 
break; 
5 
X(jJ = X_origlj]; 
if (expanded) 
deat+; 


} 


/* compute the clustering factor */ 
cf = (dea * R_1) + ea; 
if (cf < curser) { 
cur_val = value[EVAL] ; 
Eur Gr = Gis 
for (i=0; i < nvar; i++) save_coord[li] = X[il]; 


Ip 
} 
else { 

Grn=. 0; 

dea = QO; 

ea = QO; 

for (j=0; j < nvar; j++) X_origljJ = XC[jl; 
/* for each variable (direction)... */ 

for (j=0; j < nvarsiee or 1 


expanded = 0; 
/* If not on a left hand edge, move left */ 
while (X[jJ > 0) { 
6 Pi ]oss 
vcopy(vali,eval(E,X)); 
if (vali[{EVAL] && (vali[EVAL] <= value[EVAL] 
|| value[HLV])) { 
expanded = 1; 
eatt; 


else 
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break; 
} 
X({j] = X_orig(jl; 
if (expanded) { 
expanded = 0; 
deatt; 
} 
/* if we didn’t start on a right hand edge, move right */ 
while (X{jJ < R.1) f 
ee 
vcopy(val2,eval(E,X)); 
if (val2(EVAL] && (val2({EVALJ <= valueLEVAL] 
1} value({HLV])) { 
expanded = 1; 
eatt+; 
I: 
else 
break; 
} 
Kil eamkeor ie (i); 
if (expanded) 
deat++; 


} 


/* compute the clustering factor */ 
cf = (dea * R_1) + ea; 


if (!(value(HLV] && cf > TRUN)) { 
Teer serene Choe + 
cur_val = value(EVAL]; 
CursCk =e; 
for (i=0; i < nvar; i++) save_coord[i] = X{i); 


} 
} 
+} /* end if */ 
k = QO; 
} /* end while */ 


if (!last && (term == (E_orig.nterm - 1)) && !Not_all) { 
All_trun = 1; 
cur_CF = MAX_INT; 
term = -1; 
last = 1; 
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: 

if (cur_CF == MAX_INT) 
return(NULL) ; 
save_coord[nvarti] = cur_CF; 
save_coord[nvar] = cur_val; 


return(save_coord) ; 


static int valid_implicant (I) 
Implicant *1; 
[Rn nn eee eee 
SFune Elon: 
- Decide upon the validity of implicant I 
= Lecal to nace 
-clobals; 
E_work 
B_orig 
:Side_effects: 
STAT 
sCalledeby.: 
pick_implicant() 
Seales: 
next_coord() 
eval() 
vcopy () 
:returns: 
1 if a valid implicant 
O if not 


*X, 

init = i, 

Ri tv] radix =), 

value = I[->coeff, 

Vo[2] ,Vwi2]; 
static int 

coord [MAX_VAR+2] ; 


#ifdef ANALYZER 


STAT->calls_valid_implicant++; 
#endif 
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while ((X = next_coord(coord,I,init)) '= NULL) { 
eG = 10: 
vcopy (Vw, eval (&E_work,X)); 
vcopy(Vo,eval (#E_orig,X)); 
if (((VwLEVAL] < value) && !VwEHLV]) && (Vo[EVAL] < R_1)) 
return(0); 
, 


moturn(i) ; 


Static int compute_nrc(T) 
Implicant xT; 
[Rm rrr rrr erence re--- 
Sadrction: 
- Compute the NRC for the given implicant 
= Local to nd.c 
iglobals: 
radix 
nvar 
:side_effects: 
SILAT 
:called_by: 
pick_implicant() 
-calls: 
next_coord() 
eval() 
vcopy() 
sEOturns: 
- an integer NRC 


IG e *X; 
ont I_value = I->coeff; 
register 1; 
int value[2], 
Reale = Gadix — 1, 
neighbor_value[2], 
good, 
bad, 
ditt 
equal, 
neig_boun, 
fins. 
NGGa= 0. 


39 


init. = oie 
static int 
coord [MAX_VAR+2] ; 


#ifdef ANALYZER 
STAT->calls_compute_nrct+; 


#endif 
/* for each coordinate in the implicant ... */ 
while ((X = next_coord(coord,I,init)) != NULL) { 
init = ©; 
equal = 0; 


vcopy(value,eval(&E_work,X)); 
if (value({EVAL] == radix) 


continue; 
diff = value(EVAL] - I_value; 
firse = 15 
/* for each direction ... */ 
for (i=0; i < nvar; it+) { 
good = QO; 
bad = QO; 
if C(diff <=70) 928 first) 
good = 2; 
first = 0; 
} 
/* if there is a left neighbor, examine it */ 
if (x{i] != 0 &&@ x{i] == I->B[i].lower) { 
Ni aL] ase 


vcopy (neighbor_value, eval (ZE_work ,X)); 
neig_boun = neighbor_value(EVAL] - value([EVAL] ; 
X(Cil++; 
if (neighbor_value[EVAL] != 0) { 
if (!neighbor_value(HLVJ] || !value(HLV]) f{ 
if (neighbor_value[EVAL] < diff) { 
if (neighbor_value[HLV] ) 
good += 1; 
else 
bad += 2; 
sr 
if (neighbor_value[EVAL] > diff) { 
if (!neig_boun) 


bad += 2; 
if (neighbor_value(HLVj && neig_boun < 0) 
bad += 2; 


if (diff > 0 && neig_boun) { 
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if (value([HLV]) 


good += 1; 
else 
bad += 2; 
i 
} 
else { 
if (neighbor_value[HLV] || value[HLV]) 
good += 1; 
else 
good += 2; 
} 


} 
e/eend 1f */ 
} 


/* if there is a right neighbor, examine it */ 
if (XCi] '= R_1 && X({iJ == I->BLi].upper) { 
XCil++; 
vcopy (neighbor_value,eval (&E_work,X)) ; 
neig_boun = neighbor_value[EVAL] - value[EVAL] ; 
KialSs 
if (neighbor_value[EVAL] != 0) { 
if (!neighbor_value[HLV] [| !value[HLv]) { 
if (neighbor_value[EVAL] < diff) { 
if (neighbor_value[HLV] ) 
good += 1; 
else 
bad) += 2; 
} 
if (neighbor_value[EVAL] > diff) { 
if (!neig_boun) 


bad += 2; 
if (neighbor_value[HLV] && neig_boun < 0) 
bad += 2; 


if (diff > 0 && neig_boun) { 
if (value[HLV]) 


good += 1; 
else 
bad += 2; 
I; 
} 
else { 
if (neighbor_value(HLV] || value(HLV]) 
good += 1; 
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else 
good += 2; 


i 


/* update the nrc */ 
nrc = (nrc - good) + bad; 
} /* end for */ 
} /* end while */ 
return iso. 


} 
Static Implicant *pick_implicant(X) 
ane *X; 
OR mr re eer 
ounce lon: 
- Pick the best implicant for minterm X 
:globals: 

radix 
:side_effects: 

STAT 
:called_by: 

N_D() 
~Galls. 

init_implicant() 

gen_bounds () 

next_implicant() 

eval () 

vcopy() 

compue_nrc() 

copy_implicant() 

valid_implicant() 
sreturns:. 

- A pointer to a term representing the best implicant. 
a * / 
{ 
int cur_nrce = MAX_INT, 

nre = QO, 

I_value, 

i, 

init = 15 

first = 1; 


Implicant 

*I; 

Stacie mnt 

coord [MAX_VAR+2] ; 
static Bound 

I_bound [MAX_VAR+2] ; 
Static Implicant 


I_best ; 
Bound *B; 
rat vei) 
value(2]; 


#ifdef ANALYZER 
SiAt—> calls pick _implicant++; 
#tendif 


I_best.B = I_bound; 
init_implicant(X) ; 
B = gen_bounds(X) ; 
vcopy(V,eval (&E_orig,X)); 
while ((I = next_implicant(B)) != NULL) f{ 
if (V[HLV]) { 
for (I->coeff=X(nvar]; I->coeff < radix; (I->coeff)++) f{ 
if (valid_implicant(I)) { 
nre = compute_nrc(I); 
if -(farse) 
nee = 2; 
else 
Hee += 92: 
fi (nrc e<= cur nre) + 
GubencsGa- NnGEes 
t-2nrG = nrc; 
copy_implicant(&I_best,I); 


uw 


} 
first = 0; 
i 
else { 
I->coeff = X[nvar]; 
if (valid_implicant(I)) f{ 
nrc = compute_nrc(I); 
Piectirst) +. 
first = 0; 
be WGggatel << 1). ), 
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else 
Nre.+— 2: 
if (nee. <= cur_nre) 4 
cur_nre = nrc; 
l= ne = tre. 
copy_implicant(&I_best ,I) ; 


} 


returnce) past). 


dt 


APPENDIX B: AVERAGE NUMBER OF 
OUTPUT PRODUCT TERMS 


TABLE B.1: Two-Variable Four-Valued Average Output Product Terms 
Number of | Pomper and Armstrong | Dueck and Miller Neighborhood 
pat ee eee [osc 

2.109 
35575 1 
3.675 
3.934 
Seti 
3.950 
3.872 
3.699 
3.034 
3.398 
30k 
3.042 
Zeon t 
2.658 
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TABLE B.2: Two-Variable Five- Valued Average Output Product Terms 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
Input Terms Decoupling 
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TABLE B.3: Three-Variable Four-Valued Average Output Product Terms 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
=. Terms -. 





TABLE B.4: Four-Variable Four-Valued Average Output Product Terms 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
—- Terms ae 
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APPENDIX C: AVERAGE COMPUTATION 
TIME FOR MVL MINIMIZATION 


TABLE C.1: Two-Variable Four-Valued Average Computation Time 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
——F Terms Decoupling 


0.03309 
0.05517 
0.08384 
a ak o7 
0.14140 
0.16920 
0.19713 
0.22143 
0.24440 
0.27069 
0.28888 
0.31570 
0.33116 
0.34269 
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0.06156 
0.10724 
0.16931 
0.23128 
0.29612 
0.37093 
0.43512 
0.48732 
0.55478 
0.61506 
0.66258 
0.71373 
0.76289 
0.79065 


0.06336 
Cale 21 
0.17947 
0.24318 
0.31576 
0.38907 
0.46649 
0.52141 
0.59490 
0.65662 
lst 
0.77805 
0.82645 
0.85818 





0.05604 
0.09554 
0.15442 
0.22760 
0.29088 
0.36180 
0.42612 
0.50258 
0.57016 
0.62940 
0.70374 
0.76936 
0.84158 
0.88608 
0.93786 
0.99650 
1.03754 
1.07786 
1.13426 
Waar ers: 
1.18984 
e224 
1.24788 


50 


0.10954 
0.20432 
0.33730 
0.50900 
0.67988 
0.86232 
1.08536 
1.24780 
1.44212 
1.64508 
1.83080 
2.07240 
2.23954 
2.36738 
2.90298 
2.63846 
2.72054 
2.84114 
2.94840 
7 eee 
3.04716 
3.10796 
3.24612 


TABLE C.2: Two-Variable Five-Valued Average Computation Time 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
es Terms Decoupling 


OLIN 
0.20372 
0.33034 
0.50168 
0.66278 
0.84642 
1.05530 
1.23926 
1.41004 
1.60892 
1.79734 
1.98956 
2.14072 
2-24 0 
2.45116 
2.62578 
2.14588 
2.86624 
2.98710 
3.08048 
3.10488 
3.20222 
sn lZe 





0.11096 
0.17684 
0.29150 
0.43990 
0.65294 
0.86260 
1.07620 
1.36058 
133 2 
1.97194 
2.24472 
2.61166 
3.00044 
3.40472 
3.80686 
4.22300 
4.62690 
5.10534 
5.51390 
o.01o20 
6.35928 
6.69882 
el feet 
7.60032 
7.77850 
8.32284 
8.95526 
9.04034 
9.11648 
9.48198 
9.88210 
10.1734 
10.3463 
10.8508 
11.0079 
il 13s 
Piao | 
eal 
11.8719 
12.4192 
24s) 
| bases oy) | 
12.6780 


0.25190 
0.41760 
0.67308 
1.04038 
1.49640 
2.09008 
2.61164 
3.29064 
3.97860 
4.71946 
5.93036 
6.50490 
7.96094 
8.50130 
9.56092 
10.6681 
11.9961 
12.9534 
14.0393 
15.2872 
16.1587 
ests 
18.2629 
19.5004 
20.4014 
21.5390 
22.5366 
23.3190 
24.1279 
24.2701 
25.3621 
26.4659 
26.5370 
27.1308 
28.1984 
28.6137 
ooo 124 
30.3718 
30.1055 
30.7879 
31.1774 
30.4067 
31.7510 


TABLE C.3: Three-Variable Four-Valued Average Computation Time 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
Input ame Decoupling 


0.23496 
0.39676 
0.66464 
1.02516 
1.47802 
2.07168 
2.09776 
3.24702 
3.93274 
4.73014 
5.90418 
6.37408 
vga 
§.33298 
ee ea 
10.1814 
Ie 1 20 
12.0204 
sO 
14.0559 
15.1375 
15.9774 
17.0413 
17.6663 
18.4309 
19.5850 
20.7064 
21.0932 
21.9608 
22,397 I 
23.1439 
24.2141 
24.3088 
25.6463 
25.9083 
25.8980 
2(@t013 
28.0957 
28.1734 
29.1146 
29.4338 
29.0054 
30.1296 





0.35778 
0.57068 
0.89762 
1.26432 
ooo (0 
Zo 
3.20296 
4.23814 
5.21906 
6.31822 
1.86472 
9.35454 
10.9688 
13.0702 
15.9150 
17.5832 
20.8329 
22.9053 
26.1380 
30.5844 
34.5081 
39.0754 
42.8344 
47.3552 
93.0040 
Of ole 
60.7629 
63.8871 
70.0774 
13.9980 
78.4594 
84.8676 
92.5689 


O2 


1.01774 
1.64304 
2./4510 
3.97440 
5.14632 
7.31790 
9.42768 
12.6106 
15.5529 
18.5009 
ZZ Go2 
29.1544 
29.2464 
34.6918 
39.5156 
45.6882 
Szezla 
58.4218 
65.7590 
14.2779 
82.4673 
92.3658 
104.355 
111.323 
S279 
125.230 
136.164 
146.791 
160.949 
172.364 
185.003 
200.174 
21020 


TABLE C.4: Four-Variable Four-Valued Average Computation Time 


Number of | Pomper and Armstrong | Dueck and Miller | Neighborhood 
Input Terms Decoupling 


0.79966 
1.33264 
2.14 22e 
3.24550 
4.58174 
6.54000 
8.91408 
11. 382% 
13.8771 
17.0987 
20.9875 
24.6584 
29. 01g2Z 
34.6770 
40.0190 
46.5026 
54.0485 
59.8628 
67.4621 
76.7863 
86.0243 
94.9414 
100.501 
106.223 
113.879 
119.541 
NADAS 
140.376 
154.421 
165.895 
i TEO) 3 
191.966 
207.516 
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